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The iron pnictides and the cuprates represent two families of materials, where strong antiferro- 
magnetic correlation drives three other distinct ordering tendencies: (1) superconducting pairing, 
(2) Fermi surface distortion, and (3) orbital current order. We propose that (l)-(3) and the anti- 
ferromagnetic correlation are the hallmarks of a class of strongly correlated materials to which the 
cuprates and pnictides belong. In this paper we present the results of the functional renormalization 
group studies to support the above claim. In addition, we show that as a function of the interlayer 
hopping parameter, the double layer Hubbard model nicely interpolates between the cuprate and 
the iron pnictide physics. Finally, as a check, we will present the renormalization group study of a 
ladder version of the iron pnictide, and compare the results to those of the two-dimensional model. 
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I. INTRODUCTION 

The discovery of high-temperature superconductivity 
in the iron pnictides adds a new class of compounds to 
the high-Tc family^"^. In many regards, the pnictides are 
very similar to the cuprates. For instance, the undoped 
compound shows antiferromagnetic (AFM) long-range 



order and becomes superconducting (SC) upon doping, 
and the doping-temperature phase diagrams is similar to 
the cuprates. On the other hand there are important dif- 
ferences. For example, the undoped cuprates are Mott 
insulators while the undoped pnictides are semi-metals. 

In addition to the cuprates and the pnictides, there are 
other materials that become superconducting upon exit- 
ing the antiferromagnetic phase. Examples include the 
heavy fermion compounds^ and the organic compounds^ 
(in these systems pressure replaces doping as a tuning 
parameter). 

An important question concerning these materials is 
"does antiferromagnetic correlation have anything to do 
with the superconducting pairing" . In this paper, we try 
to answer this question by performing functional renor- 
malization group (FRG) calculations. We shall focus on 
the iron pnictides and the overdoped cuprates. The rea- 
son for focusing on the overdoped cuprates is because we 
believe FRG has a better chance to succeed. As for the 
iron pnictides, the general belief is that its electronic cor- 
relation is weaker. We believe this makes FRG suited for 
the entire doping range. 

Although the pnictides is not as strongly correlated 
as the cuprates, in some sense it is more difficult. In 
the cuprates there is a charge gap, below which superex- 
change generates the AFM interaction. With the AFM 
interaction, superconducting pairing appears in mean- 
field theory^. For the iron pnictides there is no charge 
gap, and the bare (purely repulsive) Hamiltonian does 
not show any mean-field SC instability. In this case 
we need a renormalization group procedure (in contrast 
to the one-shot second order perturbation procedure for 
the cuprates) to generate the effective Hamiltonian. Of 
course, the cuprate problem is difficult in another aspect, 
namely, at low energies there is the no-double occupancy 
constraint. 

For the cuprates it is well recognized that it is the 
underdoped limit which exhibits the most unconven- 
tional behaviors. Examples include pseudaogap and the 
dichotomy between the nodal and antinodal electronic 
excitations*. In that limit the system also shows the 
propensity toward charge/spin density wave order^"^^ . In 



2 



contrast the overdoped limit appears more conventional. 
The normal state looks like an ordinary metal, and the 
SC state appears to be a garden variety of d-wave su- 
perconductor. Interestingly for the LSCO system it is 
shown that even the overdoped samples show substan- 
tial AFM correlation-'^^. In the literature, the proximity 
to the Mott insulating limit is widely regarded as why 
cuprates show high Tg. In our opinion the iron pnictides 
cast some doubt on this point of view. 

Now let us narrow our focus to the recent theoreti- 
cal developments concerning the pnictides. Prom a the- 
oretical standpoint what makes the pnictides more com- 
plex is the fact that there are five relevant bands^^"^^. 
They give rise to multiple Fermi surfaces (FS). It is con- 
jectured early on that the Umklapp scattering between 
the electron and the hole FS is responsible for the SC 
pairing, and because the sign of the Umklapp scattering 
is positive, the gap function takes opposite sign on the 
electron and hole pockets^^'^'*. This conclusion is inde- 
pendently reached by a number of other groups^^"^^ by 
performing calculations with different degree of approx- 
imation. Among them the most unbiased treatment is 
the FRG of Ref.^®. This pairing symmetry and mech- 
anism is supported by recent ARPES experiments^°'^^. 
In addition, the lack of the Hebel-Slichiter peak in NMR 
experiments^^ and the presence of a neutron resonance 
peak^'^ in the SC state are both consistent with the s± 
pairing symmetry^^. 

The above approaches start from the band semi-metal 
are sometimes regarded as the "weak coupling" approach. 
While this phrase is justified for some of the above ap- 
proaches, it is not for the FRG. Indeed, as a function 
of decreasing energy cutoff various scattering amplitudes 
grow in strength. By the time only a thin shell around 
each Fermi surface is left, those scattering can get very 
strong. In addition, the cjuasiparticles contained in these 
thin shells are no longer bare electrons. For example the 
coupling of these quasiparticles to the staggered magnetic 
field is enhanced compared to bare electrons. 

There is also another approach which starts with a 
Ji — J2 model which can produce the observed (tt, 0) or 
(0, tt) AFM ordering and the associated four-fold rota- 
tion symmetry breaking^'''^^^. Upon doping, mean-field 
theory similar to those of the cuprates (except there is 
no occupation constraint) gives rise to the same out-of- 
phase SC pairing^^. However, since for iron pnictides 
there is no charge gap to justify the superexchange, the 
validity of the effective model requires a justification. In 
addition to the above, there are many other theoretical 
works addressing various aspects of the pnictides^^. 

For the cuprates and pnictides antiferromagnetism and 
superconductivity are incontrovertible. Moreover, strong 
experimental evidences indicate that the propensity to- 
ward unidirectional charge/spin density wave ordering 
(or "stripe ordering") exists for many underdoped (or 
sometimes even overdoped) cuprate families^^. Theoret- 
ically the staggered orbital current ordcr^'^° and the ten- 
dency toward Pomeranchuk instability (PI)'^'' has also 



been discussed for the cuprate materials. For example, 
the "D-density wave" order has been proposed to gen- 
erate the anti- nodal pseudogap^. Moreover, there are 
proposals that such type of order will be exhibited in the 
vortex core of the high Tc superconducting state. -'°'-'-^. 
Conc;c!rning the PI it has been proposed that interaction 
which favors a d-wave FS distortion exists^^. Moreover 
it is argued that when such distortion occur, it will lead 
to an electronic nematic state^^ which is the precursor 
of the stripe ordering. In this paper we shall show that 
a similar group of ordering tendencies exist for the pnic- 
tides. We propose that the quartet of ordering tendencies 
including the AFM ordering, the SC pairing, the orbital 
current ordering, and the Pomeranchuk instability, are 
the hallmark of a whole class of correlated materials - 
materials where the short range electronic repulsive in- 
teraction induced strong AFM fluctuations dominate the 
low temperature properties. 

The paper is organized as follows. In Sec. II, we shall 
start with the FRG results for the iron pnictides. The 
goal is to identify the interactions that grow at low ener- 
gies, and to construct an effective Hamiltonian describ- 
ing them. We shall argue that the AFM scattering is 
what's driving these strong interactions. In Sec III, we 
shall revisit the FRG result for the single-band Hubbard 
model proposed for the cuprates^^''^^. Same efforts will 
be spent to identify growing interaction at low energies. 
Many results in this section have already been obtained 
before'''' We present them for the purpose of mak- 
ing comparison with the pncitides results. It is from this 
comparison a coherent picture emerges. In this picture 
AFM scattering is the root of SC, PI, and orbital cur- 
rent ordering tendencies. In Sec IV, we show that as 
a function of the inter- layer hopping parameter (t^) the 
double-layer Hubbard model nicely interpolates between 
the physics of the pnictides and cuprates. In particular, 
as a function of increasing the band structure of this 
model evolves from two large FS (with shape similar to 
that of the curpates) to two separate electron and hole 
pockets (similar to the pnictides). Upon doping the SC 
pairing undergoes a phase transition from d-wave to out- 
of-phasc s-wavc as a function of t^. In Sec V, we present 
the renormalization group calculation for a ladder version 
of the pnictide model. The purpose of this calculation is 
to provide a check of the two-dimensional FRG results 
presented in Sec II. This is motivated by the fact that 
in the case of the cuprates, much of the d-wave pairing 
physics is already captured by the two- leg ladder^^. Sec- 
tion VI is the summary and conclusion. 

II. THE FRG RESULTS FOR THE IRON 
PNICTIDES 

The goal of the FRG calculation is to generate the 
renormalized electron-electron scattering Hamiltonian at 

low energies: 

l/(ki,k2,k3,k4,A)i/)^^_^VL,s'V'k3,s'V'k4,s- (1) 



3 




FIG. 1: (color on-line)(a) The LDA FS of iron pnictides^^ . 
1 - 3 are hole FS and 4, 5 label the electron FS. (b) The LDA 
dispersion along the (0, 0) — (tt, 0) line. Here the red, black, 
and blue curves represent the a, /3 and 7 band, respectively. 
The dashed green line marks the Fermi level of undoped com- 
pound. The Van Hove and the Dirac band singularities dis- 
cussed in the text are indicated. 



Here A is the running energy cutoff, and k = {kx,ky, a) 

specifies the momentum k = {kx, ky) as well as the band 
index a. We watch for those growing V(ki, k2, ks, k4. A) 
as A decreases. 

For the iron pnictide there are five relevant bands^^"^^. 
Two hole bands, denoted by a and /3, give rise to two hole 
pockets around the (0,0) point. Depending on the dop- 
ing, there can be another hole pocket around (7r,7r). (In 
this paper we use the so-called unfolded Brillouin zone.) 
The electron band, denoted by 7, gives two electron pock- 
ets: one around (tt, 0) and the other around (0, tt). For 
the band structure given by Kuroki et al}^ , the electron 
dispersion along (0, 0) — (tt, 0) and the FS are shown in 
Fig. 1. To a good approximation, the FS are nested by 
momenta (tt, 0) and (0, tt) for a range of doping. A mea- 
sure of the the "failure of nesting" can be obtained as 
follows. Take the best nested pair of electron and hole 
FS, displaces the hole FS by (tt, 0), record the energies of 
the electron band around the displaced FS, and calculate 
the maximum absolute value A,„is of these energies. For 
the band structure of Kuroki et al}^ ^ Amis < O.OSeF for 
< 10% hole doped or parent compound, and < Q.leV 
for < 10% electron doped compounds. For A > Amis we 
should regard the FS as nested. Above Amis the generic 
"on-shell" scattering processes include Cooper scattering, 
forward scattering, and (tt, 0), (0, tt) spin density wave 
(SDW) and charge density wave (CDW) scattering. 

For the processes mentioned above, we can cast the 
corresponding four-fermion scattering Hamiltonian into 
following form 

^^(k,p,A)Oj,6p (2) 
k,p 

where Ok can be pairing operator, current density oper- 
ator or spin density operator. (We shall explain this in 
detail later.) We then decompose ^(k, p,A) into differ- 



ent eigen-modes as 

F(k,p,A)-^u;,(A)/"(k)r(p) (3) 

i 

We monitor Wi as a function of the energy cutoff A. If Wi 
is negative and its absolute value grows as A decreases, it 
indicates that the system has a propensity for developing 
the order characterized by the following order parameter 

o, = 5]r(k)(Ok) (4) 

k 

where the form factor /'(k) tells us about the symmetry 
of the order parameter. Since the renormalized quasipar- 
ticles (left in thin shells around the FS) are not the bare 
electrons, we can not unambiguously determine the real 
space pattern of the order. The real space patterns that 
we shall present below are the simplest representation 
that captures the symmetry of /'(k) when it is Fourier 
transformed to momentum space. 

The Hamiltonian we shall perform FRG upon is given 

by 

5 

k.s a,b— 1 
i a 

a<b a<b,s.s' 

+ JHYi'4'ja'^'4'^al^^bl'^P^M +h.C.)y (5) 
a<b 

The parameters used in constructing KatO^) can be found 
in Ref.^'^. The results in this section is obtained by using 
Ui = AeV and U2 = 2eV and Jh = O.TeV. The results 
are qualitatively the same for other choices of bare inter- 
actions. The details of our FRG calculation can be found 
in Ref.i6. 



A. The Spin Density Wave 

In reference to Eq. (1) we consider the following two 
types of scattering processes, one is fci — Q = ^3 = k 
and k2 ~ Q — ki = and the other is fci — Q = fc4 = fc 
and ^2 — Q = fca = p, where Q can be either (tt, 0) or 
(0, tt). Once the momenta are fixed, the band index will 
be correspondingly fixed by the conservation of momen- 
tum. For example if k is on the electron (hole) band, 
k -|- Q will be on the hole (electron) band. The corre- 
sponding scattering terms are 

F(k + Q, p + Q, k, p, A)v4+Q,,V4+Q,,'^k,.'V'p.s, (6) 
and 

Vik + Q, p + Q, p, k, A)v4+Q,,V4+Q,,'^p,.'V'k,s, (7) 



We now cast their sum into SDW and CDW scattering 
as follows: 

VsDw(k, p, A)^k,Q<S'p,Q + VcDw(k, p, A)nk^Qnp^Q, (8) 
where 

Sk,Q = i^lk+Q^ss'i^s'M- (9) 



]/ / / aI 



and 



nk,Q =Xl''^Ik+QV's,k• 



Using the identity 

^ ss' ^rr' — sr' s' ^ ss' ^tt' 

one can show that 



(10) 



(11) 



FsDw(k, p, A) = -iy(k + Q, p + Q, k, p, A) (12) 

and 

VcBwiK p, A) = y(k + Q, p + Q, p, k, A) 

-iF(k + Q,p + Q,k,p,A) (13) 

In the above and later section, by "CDW" we mean the 
(tt, 0) or (0, tt) density wave in the spin singlet channel. 
According to this definition an orbital current density 
wave is categorized as "CDW" . In the rest of this sub- 
section we shall first focus on the SDW scattering only. 

The SDW instability already exists at the mean field 
level in the bare Hamiltonian level. For the type of in- 
teraction considered in Eq. (5) the representative pattern 
of the mean-field magnetic order is shown as the type-a 
SDW in Fig. (2). It has been shown in Ref.^^ that due 
to the point group symmetry of the band wave function, 
the off-diagonal SDW matrix element between the larger 
hole pocket (labeled by 2 in Fig. (1)) and electron pocket 
(labeled by 4 and 5) vanishes on the principal axes of the 
square lattice. As a result the form factor has a d-wave 
symmetry and the SDW state is always gapless. This 
feature is reproduced by the FRG. More specifically, the 
leading channel SDW ,i.e., the channel associated with 
the most negative Wi , has a d-wave form the form factor 
as shown in Fig. 3(a). Note here, since the SDW order 
parameter involves two Fermi surfaces, each form factor 
in Fig. (3) (a) is labeled by a pair of FS indices. For a 
range of doping, this type of SDW is the dominant order- 
ing tendency at low energies. In addition, FRG shows a 
subdominant SDW order. The real space representative 
is shown as the type-b SDW in Fig. (2). In this pattern 
the magnet moment is mainly associated with the bonds. 



B. The Superconducting Pairing and its 
Mechanism 



Here we focus on ki = — k2 = k and ka = — k4 = p in 
reference to Eq. (1). These process can be expressed as 
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FIG. 2: (color on-line) A schematic real space representation 
of the two types of SDW order. 



the sum of singlet and triplet pair scattering 

ysc,s(k, p, A)A^A,,p -I- Vsct (k, p, A)At ,^At,p. (14) 
Here 

l^sc,s(k,p,A) 

= l/(k, -k, p, -p, A) + l^(k, -k, -p, p, A) 
+l/(-k, k, p, -p. A) + y (-k, k, -p, p. A) (15) 

ysc,t(k,p. A) 

= y (k, -k, p, -p. A) - F(k, -k, -p, p. A) 

-y(-k, k, p, -p. A) + y(-k, k, -p, p. A), (16) 



and 



^^■■^ " ^/2 



7f 

^(<kV'S,_k-<_k^.\k)- (17) 



V2 

For all the parameter we have looked at the triplet pair- 
ing tendency is always much weaker than that of the sin- 
glet. Among the singlet pairing, we find that the leading 
mode has s± symmetry, i.e. the pairing order parameter 
does not change sign around each FS, but takes oppo- 
site sign between hole band (a and /3) and electron band 
(7). The associated form factor is shown in Fig. 3(b)^^. 
Unlike the SDW, each FS has one form factor. 

One big advantage of FRG is that it allows one to mon- 
itor the growth of different classes of scattering processes 
as a function of the energy cutoff, so it is possible to 
identify which class grows strong first and hence poten- 
tially serves as the driving mechanism for other classes 
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FIG. 3: (color on-line) (a) The FRG form factors of the SDW 
order parameter. The red , black and blue curves are the form 
factors associated with the (2, 4), (1, 4) and (3, 5) Fermi pock- 
ets (i.e. when (k, k-|-Q) £ (2, 4), (1, 4), (3, 5) FS), respectively, 
(b) The FRG form factors of the SC order parameter. Here 
different colored curves represent the order parameters asso- 
ciated with the five different FS (i.e. when (k, — k) G 1, 5 
FS). For each pocket 6 is defined as the angle sustained by the 
vector connecting the center of the pocket to point {kx, ky) on 
the FS and the vector (7r,0). The k-coordinate of the center 
of FS are (0,0) for pocket 1,2, (tt, tt) for pocket 3, (0, — tt) 
for pocket 4 and (7r,0) for pocket 5). The parameters under 
which the above results are obtained are specified in the text. 



of scattering. As discussed in Ref.^^ we always find the 
SDW scattering (this includes a whole class of processes) 
grows strong the first. This is related to the fact that 
even the bare Hamiltonian shows a SDW mean-field in- 
stability. What subsequently triggers the SC scattering 
is a restricted set of scattering processes which have dual 
character of SDW and SC scattering as shown in Fig. 
4(a) In terms of the scattering Hamiltonian these spe- 
cial type of scattering processes are 

V{k + Q, -k + Q, k, -k, A)^At^Q_^^^li^^Q_^,V'k,.'^-k,s. 

(18) 

The dual character of this type of scattering can be seen 
as the possibility of "factor" it into both the SDW and SC 
channels. Together with other SDW scattering, this type 
of interaction grows positively large in early stages of the 
FRG (Fig. (4)(b))i^ Upon SC factorization this scat- 
tering give rise to (positive number)x A]^ , gAk, hence 




FIG. 4: (color on-line) (a) A schematic representation of the 
Umklapp scattering processes which has the dual character of 
SDW and SC. (b) The amplitude of the the scattering process 
shown in (a) as a function of the FRG energy cutoff A. (c) 
A schematic real space representation of cos kx cos ky singlet 
pairing. 



favors a SC order parameter satisfies 

(Ak)(AUQ><0 (19) 

Once the above special type of Cooper scattering is 
seeded, other members of the inter-pocket Cooper scat- 
tering processes grow strongly positive in subsequent 
FRG steps. Two real space representatives of the pairing 
satisfying Eq. (19) is are the two types of second neigh- 
bor pairing (Fig. 4(c)) with form factor coskx cos ky and 
sin kx sin ky. (Of course there are higher angular momen- 
tum pairing that satisfies Eq. 19, however they in general 
produce more nodes on the FS and hence are less favor- 
able.) The former has s± symmetry and the latter has 
dxy symmetry. However, given the locations of the FS 
of the iron pnictide, the latter is less favorable because 
it produces four nodes on each FS while the former does 
not. It is interestingly that precisely these two types of 
pairing symmetry are observed as the leading (s±) and 
one of less divergent (dxy) SC pairing instabihties^^. 

The second neighbor singlet pairing in Fig. 4(c) is 
also consistent with the short-range magnetic correla- 
tion shown in Fig. 2(a). Indeed according to Fig. (2) 
the next-nearest pairs always preferred to be anti- 
ferromagnetically correlated. This is in harmony with 
the singlet pairing depicted in Fig. (4)(c). From our cal- 
culation, the orbital content of the Cooper pair is domi- 
nated by dxz(YZ) - dxz{YZ), dxz(YZ) - dyz^xz), and 
dxY — dxY- 

The SC gap function determined from our FRG anal- 
ysis often has large variation in the electron pocket. For 
instance, in Fig. 3(b) the gap nearly touches zero along 
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one of the principle axes ((0, 0) — (tt, 0)) around the elec- 
tron FS. Such large gap variation is shown to give rise 
to power law-like NMR relaxation rate so long as the 
temperature is not smaller than the minimum gap'^''. It 
can also be consistent with the thermal transport exper- 
iments results obtained recently'^^. The origin of such 
gap variation is the variation of the orbital content of 
the band wave function around the electron pocket^®. 
According to band structure^^ the electron band eigen- 
function along (0, 0) — (tt, 0) is primarily XY-\ike. In 
contrast along (tt, 0) — (tt, tt) direction the electron band 
eigenfunction is dominated by the XZ and YZ orbitals. 
So if the orbital content of the pairing function is such 
that XZ{YZ) have more weight than the XY orbital, the 
type of gap variation will result. Moreover this would also 
predict that the gap associated with the two hole pockets 
around (0, 0) (labeled by 1 and 2) are larger that that of 
the (tt, tt) hole pocket (labeled by 3) because the former 
are mainly made oi XZ and YZ orbitals and the latter 
is mostly XY-\ike. This is indeed consistent with the 
FRG result. The above interpretation of the gap varia- 
tion is consistent with a recent FRG calculation where 
the orbital dependence of the interaction matrix element 
is ignored^". Their resulting gap function is far more 
uniform than ours. 

If by tuning some bare interaction parameters, one can 
reverse the weight of pairing function on different orbitals 
(i.e. to make the XY orbital weight more than the XZ 
and YZ) it is possible to reverse the gap anisotropy dis- 
cussed above. In addition, we also suspect that, on the 
sample surface where the four-fold rotation symmetry is 
broken, the orbitals content of the band wavefunction 
will be modified. This could give rise to a different (less) 
degree of gap anisotropy. If so this could reconcile why 
surface sensitive measurement such as ARPES observed 
a far more isotropic gap^^. 

C. The Pomeranchuk Instability 

In this section we focus on ki = k4 = k, k2 = ka = p 
or ki = ka = k, k2 = k4 = p in reference to Eq. 1. 
Similar to Eq. 12 and Eq. 13 we can decompose these 
forward scattering into the singlet and triplet channels 
as 

Fpi,,(k, p, A) = ^^(k, p, p, k, A) - ^V{K p, k, p. A) (20) 
and 

ypi,t (k, p, A) - -^Vik, p, k, p, A), (21) 

These channel decoupling allows the forward scattering 
to be written in the following Fermi liquid form 

T/pi,s(k, p, A)nknp + Fpi,t(k, p, A)S^k5p (22) 

where 

= ^ V'lk'/'sk (23) 

s 




S 3£ 2;r 1.2 0.6 0.3 0.15 0.07 



FIG. 5: (color on-line) (a) A schematic representation of the 
scattering processes which have the dual character of SDW 
and PI. (b) The amplitude of the the scattering process shown 
in (a) as a function of the FRG energy cutoff A. (c) The form 
factor, /k, in Eq. (25) associated with the most negative Wi. 
The angle 9 is defined the same way as that in Fig. 3 



and 

'5^k = 5I^IkO's.'^«'.k (24) 

s 

As before we can decompose Vpi,s,t(k, p, A) into eigen 
modes as 

yp,.,t(k,p, A) = ^m.(A)r(k)r(p). (25) 

i 

When none of the Wi are negative the FS is stable. 
However, when certain Wi become sufficiently negative 
(so that it can overcome the stiffness imposed by the 
local Fermi velocity) FS deformation occurs. Attrac- 
tion in Vpi.s will drive spin-independent FS deformation 
while attraction in Vpi.t drives spin-depedent FS defor- 
mation. In our FRG calculation, we find Vpi^t is usually 
weaker than Vpi^s- As a result, we shall focus on the 
spin-independent PI in the rest of the subsection. 

According to our FRG result, the strongest forward 
scattering occurs between quasiparticles on the electron 
and hole pockets. Analogous to the inter-pocket Cooper 
scattering, there is a special set of such forward scatter- 
ing processes that have the dual character of being SDW 
scattering, namely, 

V{k -f Q, k, k, k + Q, A)i'l+Q^,i'L^^ks'^k+Q., (26) 

as shown in Fig. 5(a). Together with other SDW scatter- 
ing this special type of forward scattering grows strongly 
positive upon FRG, as shown in Fig. 5(b). Viewed as for- 
ward scattering this scattering favors the FS deformation 
at k and k -f Q to have opposite sign (due the positive 
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FIG. 6: (color on-line) (a) The FS distortion triggered by 
the PI. After the distortion each electron pockets is split into 
two smaller electron pockets, (b) For very strong forward 
scattering each electron pockets turns into two hole pockets 
after the FS distortion. In addition the smaller hole pocket 
centered around (tt, tt) is gone. 



scattering amplitude) . Physically this amounts to the 
transfer of quasiparticles from k to k + Q or vice versa. 
In terms of Eq 25 this amounts to 

/(k)/(k + Q) <0. (27) 

Similar to the "seeding effect" discussed in the last sub- 
section, this special set of forward scattering triggers 
other forward scattering processes to grow in subsequent 
renormalization group steps. According to our FRG 
analysis the /'(k) associated with the most negative Wi 
also has s± symmetry (see Fig. 5(c)), consistent with the 
requirement of Eq. (27). 

A strong s±-channel forward scattering favors the 
transfer of electron from the electron to the hole pock- 
ets or vice versa. As a result, it causes the electron and 
hole pockets to either both shrink or both expand. The 
band curvature further determines which possibility is 
energetically more favorable. As the pockets expand, the 
Fermi velocity gets larger, hence the kinetic energy cost 
increases with the deformation amplitude. The reverse 
is true for the shrinking Fermi pockets. As the result the 
band curvature always favor both electron and hole pock- 
ets to shrink. For a semi-metal with an even number of 
electrons per unit cell, and simple quadratic dispersing 
bands, strong repulsive inter-pocket forward scattering 
will trigger a transition to the band insulator. 

The situation for the iron pnictides is more complex 
due to the band structure. Shghtly beneath each of the 
electron FS there is a Van Hove singularity. A little bit 
below that there are two Dirac points (two for each elec- 




FIG. 7: (color on-line)(a) The form factor, /k, in Eq. (25) as- 
sociated with the sub-leading Wi of Pomeranchuk instability. 
The angle 9 is defined the same way as that in Fig. 3; (b) 
The FS distortion triggered by the d^2_y2-type PI. After the 
distortion one of the electron pockets is better nested with 
the hole pocket. 



tron pocket) . In the presence of these band singularities 
the PI discussed above can induce a topological change in 
the shape of the electron pockets. As the Wi associated 
with the s± forward scattering channel gets more neg- 
ative, a first order transition where the "Fermi energy" 
(the band energy separating filled and empty states after 
the distortion) in the 7 band jumps to the vicinity of the 
Van Hove singularity. However, this does not mean a dif- 
ferent chemical potential for the electron and hole bands, 
because the fermi liquid interaction energy also needs to 
be taken into account. For still more negative Ws± all 
pockets monotonically shrink as a function of decreasing 
Ws±- Before the 7-band (the blue curve in Fig. 6(al)) 
completely empties out, two electron like pockets will ap- 
pear near each M point as shown in Fig. 6(a2). For the 
parent compound strongly enough inter-pocket forward 
scattering can empty out all the pockets and hence induce 
a transition to the band insulator. For the hole doped 
case, when the 7-band becomes completely empty, there 
are still two sizable hole pockets near the (0, 0) point. As 
the forward scattering gets even stronger, the portion of 
the /?-band (black curve in Fig. 6(bl)) beneath the Dirac 
points can be depopulated. When that happens two hole- 
like Fermi pockets will appear near each (7r,0) or (0,7r) 
point as shown in the left of Fig. 6(b2). The possibility of 
an interaction-driven modification of the electron pocket 
topology has also been studied recently in Ref.''^ where a 
three-band model and the "FLEX" approximation were 
employed. 



D. A Civ Breaking Pomeranchuk Distortion and its 
Interaction with the AFM Order 

In addition to the dominant s± distortion channel, our 
FRG also predicts a sub- leading channel that has dx2-y2 
symmetry, as shown in Fig. 7(a). This is a band version 
of the orbital order. The distorted FS breaks the crystal 
C41, symmetry. For example, since the deformation of the 
two electron pockets have opposite sign, one pocket will 
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(c) 



FIG. 8: (color on-line) (a) A schematic representation of the 
scattering processes which have the dual character of SDW 
and CDW. (b) The amplitude of the the scattering process 
shown in (a) as a function of the FRG energy cutoff A. (c) The 
form factor /k associated with the leading CDW instability. 
The angle 6 is defined the same way as that in Fig. 3; (d) A 
schematic real space representation of orbital current density 
wave order. 



expand and the other will shrink. As for the hole pock- 
ets, it wiW expand in one direction and shrink in the or- 
thogonal direction, which makes the hole pocket slightly 
elongated. When this happens, one out of the two elec- 
tron pockets will be better nested with the hole pocket, 
as shown in Fig. 7(b). As a result a choice among the two 
possible SDW ordering wavevectors ((0,7r) and (7r,0)) 
will be made, and an accompanied lattice distortion is 
expected. In fact, tetragonal to orthorhombic structure 
transition has been observed in both "1111" and "122" 
compounds^^, and theoretical explanation based on J1-J2 
model has been proposed^^. Thus there is good reason to 
believe that the C^y breaking distortion will set in as the 
SDW order develops. Recent inelastic neutron scatter- 
ing experiment has revealed a discrepancy between the 
spin wave dispersion near (tt, tt) and the prediction of the 
Ji — J2 model^^. In particular, while the Ji — J2 model 
predicts a minimum, the actual spin wave dispersion ex- 
hibits a maximum at (tt, tt). This discrepancy reveals 
the fact that the 90— degree rotation symmetry break- 
ing in the magnetic state has a substantial effect on the 
electronic structure which in turn modify the magnetic 
excitation spectrum. In our picture once the SDW-driven 
Civ breaking FS distortion occurs, the electronic struc- 
ture become asymmetric with respect to 90— degree ro- 
tation. This distortion lifts the degeneracy between the 
(tt, 0) and (0, tt) orders. So long as the asymmetry is suf- 
ficiently big we would expect the minimum in the spin 
wave dispersion at (tt, tt) to turn into a maximum. 



E. The Orbital Current Density Wave (the 
Imaginary part of CDW) Order 

As stressed in Sec. II (B) by "charge density wave" we 
mean (tt, 0) or (0, tt) singlet density wave. This include 
both the conventional charge density wave with charge 
density modulation or orbital current density wave where 
no charge modulation is present. 

There is no CDW instability in the bare Hamiltonian. 
Upon FRG, negative WiS associated with the CDW chan- 
nels are generated. Interestingly the leading CDW insta- 
bility is always of the orbital current type. Similar to the 
Cooper and forward scattering, CDW also share a set of 
dual scattering processes with the SDW: 

T/(k, k, k + Q, k -f Q, A)^i+ciji+ci,,^P^.'^Pw.. (28) 

as schematically shown in Fig. 8(a). Like all dual char- 
acter scattering amplitudes this type of scattering also 
become strongly positive upon FRG (Fig. 8(b)). This 
positive scattering favors 

dk4+Q < 0. (29) 

where 

= X!^k+QsV'ks- (30) 

Replace k by k -I- Q in Eq. (29) it is simple to show 

that rfk+Q = dk- Combining this result with Eq. (29) we 
conclude that 

dkdk < 0, (31) 

as the result an imaginary CDW order parameter (dk) is 
favorable. This implies the CDW is of the orbital current 
type. The form factor of the most divergent CDW chan- 
nel is shown in Fig. 8(c) and a real space representative 
of this type of CDW order is illustrated in Fig. 8(d). 
Similar to the SDW case, the relation between the mo- 
mentum space form factor and the real space pattern is 
not transparent. This is due to the complex momentum 
dependance of the Bloch wavefunctions. The possibility 
of CDW order has also been recently noted in Ref [^^]. 

F. Competing Order - the Grand FRG Flow 

In Fig. 9 we plot the Wi{A) (in reference to Eq. (2) and 
Eq. (3)) for all the leading scattering channels discussed 
in the last four subsections. We separate the flow into the 
high, intermediate and low energy regimes. In the high 
energy regime (i.e. A ~ bandwidth, e.g., A > LOe^), the 
effective Hamiltonian is very close to the bare Hamilto- 
nian in Eq. (5). Here the only instability is in the a-type 
SDW channel (the blue line in Fig. 9(a) and(b)). In the 
intermediate energy range, (e.g., 0.15eV < A < l.OeV), 
a number of different ordering tendencies are in compe- 
tition with one another. As shown in Fig. 9(a2) and 
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(a) SC sample (b) SDW sample 

FIG. 9: (color on-line) The grand FRG flow for the iron pnic- 
tides. Here the leading Wi{A) for the SDW,SC,PI and CDW 
channels are plotted as a function of A. (a) is constructed 
using the parameters describing a SC sample and (b) is the 
corresponding plot for a SDW sample. (al),(bl) cover the full 
energy scale and (a2),(b2) zoom-in at the intermediate energy 
regime, with the competing order manifested. 



(b2), this includes different types of SDW, SC, PI and 
CDW discussed in the last four subsections. In this en- 
ergy range their relative Wi are dependent on the ini- 
tial conditions, such as the bare interaction and doping. 
However, the fact that they are all present at the in- 
termediate energy is very robust. In the lowest energy 
range, (e.g. A < O.lbeV), the a- type SDW (blue Hne) 
and the s± SC pairing (red line) dominate over all oth- 
ers. Depends on the doping we show two cases where 
the SC pairing is the most "divergent" scattering in Fig. 
9(a) and the SDW scattering is the most divergent in 
Fig. (9)(b). The form factors for these two dominant 
order tendencies have been shown in Fig. 3. 

Thus like the cuprates (see Sec. Ill) iron pnictides 
exhibit competing order at intermediate energy scales. 
More specifically both system exhibit the same set of or- 
dering tendencies, i.e., AFM, SC, PI, and CDW. As dis- 
cussed in the introduction we view this group of compet- 
ing order as the hallmark of a class of strongly correlated 
system governed by short-range repulsive interactions. 



G. The validity of the Ji — J2 Exchange as the 
Effective Interaction for the Pnictides 

In Ref.^^, much credit to the authors, it is proposed 
that the Ji — J2 magnetic exchange interaction 

captures the AFM and SC correlation of the pnictides. In 
the above (ij) represents the nearest neighbor bonds, and 
{{ij)) represents the next nearest neighbor bonds. The 
spin operator expressed in terms of the lattice electron 
operator is given by = Y.a,s,s' clas^ss'Cias' where a are 
the three Pauli matrices and a is the orbital index. After 
Fourier transform this gives rise to the following electron- 
electron scattering Hamiltonian in the momentum space 

J2Vj,^J^ ^ki,k2,k3,k4, V'ki,sV'k2,s'V'k3,s'V'k4,s, 

(33) 

where we have factored out an overall energy scale ~ J2 . 

In this subsection we test whether such an interac- 
tion reproduces the dominant low energy interaction gen- 
erated by our FRG. The minimum requirement is for 
Eq. (33) to capture the lowest energy SDW and Cooper 
scattering. More specifically, we regard Eq. (33) as the 
fully rcnormalized interaction in Eq. (1) and extract the 
scattering in various channels discussed in subsection 
IIA-IID. We then decompose the associated scattering 
into eigen channels as in Eq. 3. We check whether (i) 
the Wi associated with the SDW and SC scattering are 
more negative than all other channels, and (ii) whether 
the symmetry of the leading SDW and SC form factor 
agrees with that of the a-type SDW and s±. If both (i) 
and (ii) are passed we compute the overlap between the 
form factor predicted by the Ji — J2 model and our FRG 
result. In practice, we use the same set of discrete k 
points on FS to determine the form factor of the Ji — J2 
interaction as used in our FRG calculation. 

For — 1.5 J2 <'--^ Ji <'--^ J2, the Ji — J2 model passes (i) 
and (ii) and the results for the overlaps are shown in Fig. 
10. For both SDW and SC the best overlap is reached 
when a small and negative Ji is allowed. Thus at the 
lowest energy, the Ji — J2 model passes with flying colors. 
For the intermediate energy range the Ji — J2 model can 
also reproduce the form factor for the spin-independent 
Pomeranchuk distortion described in subsection 11(C). 
However it can not reproduce the CDW form factor and 
other subleading SDW,SC and PI form factors. Until 
writing we have not succeeded in finding a simple local 
model that can capture the form factor associated with 
all channels. 

Thus we conclude the Ji — J2 model is sufficient to de- 
scribe the low energy physics of the iron pnictides where 
only SDW and SC are important. We do not think it 
should be used as a lattice model with the full kinetic 
bandwidth. 
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FIG. 10: (color on-line) The overlap between the form factors 
predicted by the Ji — J2 model and those obtained from the 
FRG. The black and blue curves are the overlaps associated 
with the SC SDW form factors, respectively, (a) and (b) 
corresponds to SC and SDW samples in Fig. 9. The overlap 
is set to zero when the symmetry of the form factors do not 
agree, or when SDW and SC are not the leading instability 
of the Ji — J2 model . 



III. THE FRG RESULTS FOR THE ONE-BAND 
t t' HUBBARD MODEL 

In this section we present the FRG results for i — i' 
single-band Hubbard model. This model was used in 
Ref.'^'* to describe the cuprates. It is well known that for 
the cuprates the Hubbard interaction is a bit too strong 
for the FRG approach. Here we take a different point of 
view. Let us imagine a material with exactly the same 
band structure as the cuprates but with a weaker cor- 
relation. Will this material exhibit similar ordering ten- 
dencies as the cuprates? If so, it can be argued that this 
material is adiabatically connected to the cuprates. In- 
tuitively one expects this type of intermediate coupling 
system to resemble more the overdoped than the under- 
doped cuprates. The main purpose of this section is to 
compare the FRG results for this type of hypothetical 
system with those of the iron pnictidcs discussed earlier. 
The hope is that some hints on the pairing mechanism 
for both will emerge from such a comparison. 

The Hamiltonian we studied here is defined as 
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FIG. 11: (color on-line) The FRG resuhs for the t - t' Hub- 
bard model, (a) The SDW form factor and a schematic real 
space representation of the (tt, tt) AFM order, (b) The d-wave 
SC pairing form factor and its schematic real space represne- 
tation. (c) The PI Form factor and the distorted FS surface. 
The blue dashed line and the red solid line denote the original 
and the distorted FS, respectively, (d) The imaginary CDW, 
or the orbital current order, form factor, and its schematic 
real space representation. 



W = ^ ek^ki + '^{U ni^n^i - ^xm) (34) 

k,s i 

where the band dispersion is 

ek = — 2t(cos kx + cos ky) + At' cos k^ cos ky — fi. (35) 

Hereafter in unit of t (i.e. t = 1) we shall take t' = 0.2 
and [/ = 4, we shall set /i = —0.7 which corresponds 
to 10% hole doping. The band dispersion ek of Eq. (35) 
gives a FS shown as the blue curves in left panels (a)-(c) 
of Fig. (12). 



A. D-symmetry SC Pairing, FS Distortion and 
Orbital Current Order 

Upon FRG the (tt, tt) SDW and the Cooper scattering 
grow strong at low energies. Their corresponding leading 
channel form factors (in reference to Eq. (3)) are shown in 
Fig. ll(a,b). The SC pairing factor factor has the well- 
known dx2-y2 symmetry. In addition, the FRG also gen- 
erates the PI and (imaginary) CDW ordering tendency 
at low energies. Their associated form factor also have 
the symmetry. The d^^-y^ PI will trigger a C^y- 
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-T -| I ;r A/t 

(C) SDW/CDW 

FIG. 12: (color on-line) The dual scattering processes for the 
one-band t - t' Hubbard model, (a) The dual SDW-SC pro- 
cesses, (b) The dual SDW-PI processes, and (c) The dual 
SDW-CDW processes. The left column is the schematic of 
the scattering processes and the right column shows the am- 
plitude of the the scattering process as a function of the FRG 
energy cutoff A. 



breaking FS distortion as shown in Fig. ll(c)'^-'^. The 
above CDW order is the famous "DDW" orbital current 
state (Fig. ll.(d)). 

Like the pnictides the strong SC, PI and CDW scatter- 
ing all grow out of the dual scattering processes shared 
by them and the SDW scattering . For example, the 
scattering Hamiltonian describing the dual SC and SDW 
scattering is given by Eq. (18) (here Q = (tt, tt)) and is 
illustrated in Fig. 12(a). Since such Cooper scattering is 
shared with SDW, it has a positive amplitude, and favors 
(Ak)(Ak+Q) < 0. Two representative real space pairing 
form factor whose Fourier transform has this symmetry 
are given by cos + cos ky and cos k^ — cos ky . For small 
t' /t the node lines of cos + cos ky (cos kx + cos ky — 0) 
lies very close to the FS and intersects it many times. On 
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FIG. 13: (color on-line) (a) The grand FRG flow for the 
t — t' Hubbard model. Here the Wi{A) associated with the 
SDW,SC,PI and CDW are plotted. The inset zooms in at the 
intermediate energy regime, (b) The overlap of the FRG form 
factors with those predicted by the Ji — J2 model. The red 
and the blue curves are the overlaps for the SC and SDW form 
factors, respectively. The overlap is set to zero when the sym- 
metry of form factors do not agree or when the Ji — J2 model 
does not predict SC and SDW as the leading instabilities. 



the other hand the node lines of cos k^ — cos ky basically 
runs perpendicular to the FS and only intersects it at four 
points. As a result the latter pairing symmetry is more 
energetically favorable. Similarly, the PI and SDW share 
the processes of Eq. 26 as shown in Fig. 12(b). These 
repulsive processes favor the transfer quasi-particle from 
k to k -I- Q and yield a FS deformation as shown in Fig. 
11(c). CDW and SDW share the process of Eq. (28) 
and it is shown in Fig. 12(c), which gives rise to charge 
current order for the same reason as discussed in section 
II (D). 

We show in Fig. 13 the leading w, for the SDW, SC, PI 
and CDW scattering as a function of the energy cutoff A. 
Similar to the iron pnictides case only SDW channel has 
negative Wi for large A. For intermediate A, as shown in 
the inset of Fig. 13, the Wi associated with all the scat- 
tering channels mentioned above are negative, signaling 
competing order at intermediate energies. For small A 
SC (and SDW) ordering tendencies clearly surpass all 
others. Following the same strategy as used in Sec. II, we 
fit the FRG form factor of SC and SDW to those pre- 
dicted by a Ji — J2 model. As shown in Fig. 13(b), once 
J2 becomes positive, the fit worsens immediately. While 
the fit remains very good for a large range of negative J2 
until some other instability becomes dominating. This is 
because a negative J2 does not introduce any frustration 
to the (tt, tt) SDW correlation. 

B. The Stripe Tendency 

The propensity to stripe order has been observed by 
several experimental probes (e.g., inelastic neutron scat- 
tering, and STM) for several families of underdoped 
cuprates. In particular inelastic neutron scattering has 
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FIG. 14; (color on-line) (a) A schematic illustration of the 
distorted FS and the scattering processes that have the dual 
character of SC and incommensurate uni-direction charge 
density wave, or in short stripe order. The blue solid line 
is original FS and the red dashed line is the distorted FS trig- 
gered by the PI. (b) The FRG flow of the scattering amplitude 
associated with (a). 




FIG. 15: (color on-line) (a) A schematic illustration of hop- 
ping associated with the double payer Hubbard model, (b) 
A schematic representation of the proposed phase diagram, 
where S is doping. 



IV. THE DOUBLE-LAYER HUBBARD MODEL 



found incommensurate peaks at low energies for LSCO 
and YBCO^^. At high bias (comparable to the pseu- 
dogap energy scale) STM has revealed glassy pattern of 
local 90-degree rotation symmetry breaking, and a spa- 
tial modulating local density of states with a period of 
- 4-5 lattice spacing in Bi2212 and NaCCOC'*. Recent 
ARPES experiments performed on Bi2212 have found 
the evidences for two different gaps. The gap in the 
anti-nodal region appears to be particle-hole asymmet- 
ric consistent with a density wave gap^^. These experi- 
mental evidences point to the possibility that the gap in 
the antinodal region is due to the formation of a density 
wave, while that in the nodal region is due to SC pairing. 

In this subsection we examine the possibility of stripe 
order from the FRG point of view. It has been pointed 
out early on'* that the nearly parallel (hence well nested) 
FS segments in the antinodal regime could cause the 
propensity for CDW order. The nesting wavevectors are 
(q, 0) or {0,q) with q ~ 27r/5. Interestingly, after the 
PI distortion, the antinodal FS become better nested 
(Fig. (14)) in one direction and the nesting vector is 
slightly increased. Interestingly the scattering process 
shown in Fig. 14(a) is a also Cooper process. When the 
d-wave Cooper scattering becomes very strong, the am- 
plitudes associated with Fig. (14) (a) becomes strongly 
negative (because it connects pair states with the same 
gap sign). Upon further FRG the scattering amplitudes 
associated with these SC/stripe dual processes continu- 
ously to grow in magnitudes (Fig. 14(b)). When this 
type of scattering is decoupled in the CDW channel, due 
to its negative amplitude, favors a real CDW order pa- 
rameter. This implies real charge density modulation 
rather than non-zero orbital current. The uni-directional 
charge density wave resulted from this could be consis- 
tent with charge stripes. 



As we emphasized in the previous sections, from the 
FRG point of view, the key prediction of the magnetic 
scattering driven SC is that the SC order parameter has 
to satisfy AkAk+q < 0. In general there could be many 
different pairing functions satisfying this constraint. Fi- 
nally which pairing function is more energetically favor- 
able is determined by which form factor can give rise 
to maximal gapping of the FS. In this section, we will 
further strengthen this point by studying a double-layer 
Hubbard model which is schematically shown in Fig. 15. 
It consists of Hubbard layers coupled by a vertical tun- 
neling tz'*^. The Hamiltonian is given by 

n^Hl+H2~t,Y^ C^,^,.C2,^,. + h.C. (36) 

Within each layer, it is described by the Hamiltonian 
Hi^2 as in Eq. (34)(except we set t' = for this study). 

With the interlayer hopping, the vertical pair of sites 
in each unit cell form bonding and anti-bonding combi- 
nations as follows 

d2,Ls = {C1A,S - C2,i,s) • (37) 

In terms of the d operators the Hubbard interaction be- 
comes 

i.s,s' 

i 
i 

+ "2" X] 4,2,s^l,i,s'^2,z,s'^^2,2,s + 4,^,s4,^,s'^l,^s'C^l,^,s■ 
(38) 
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FIG. 16: (color on-line) (a) The FS of the double-layer Hub- 
bard model with tz — 0.5 at 4% hole doping. The two solid 
lines (blue and red) are the electron and hole FS respectively. 
The purple dotted line is the nodal line of cos kx — cos ky , and 
the black dashed line is the nodal line of cos k^ + cos ky. (b) 
The Wi{A) associated with SDW, extended s-wave SC and 
d-wave SC. (c) The FRG predicted d-wave pairing form fac- 
tor, (d) The FRG predicted form factor of the subleading 
extended s-wave pairing. Here for electron FS around (0, 0) 
, 8 is the angle between {kx,ky) and (7r,0). For the hole FS 
around (7r,7r), 6 is the angle between (k^ — n,ky — tt) and 
(^,0). 
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FIG. 17: (color on-line) (a) The FS of the double-layer Hub- 
bard model with tz — 1.0 at 3% hole doped (b) The Wi(A) 
associated with SDW, extended s-wave SC and d-wave SC. 
(c) The FRG predicted d-wave pairing form factor, (d) The 
FRG predicted form factor of the nearly degenerate extended 
s-wave pairing. The angle is defined in the same way as 
Fig. (16). 



FIG. 18: (color on-line) (a) The FS of the double-layer Hub- 
bard model with tz = 1.5 at 8% hole doping, (b) The Wi{A) 
associated with SDW, extended s-wave SC and d-wave SC. 
(c) The FRG predicted extended s-wave pairing form factor, 
and (d) (d) The FRG predicted form factor of the less diver- 
gent d-wave pairing. The angle is defined in the same way as 
Fig. (16). 



The energy bands formed by the bonding and anti- 
bonding orbitals are given by e(k) ± tz where e(k) = 
— 2t(cos kx + cos ky) + At' cos cos ky — fj,. Thus, as long 
as tz is smaller than the band width of e(k), there will 
be an electron FS around (0, 0) and a hole FS around 
(tt, tt). These FS are well nested at half- filling (two elec- 
trons per unit cell) since their size are the same. Thus 
with any non-zero positive U we expect Q ~ (tt, tt) 
SDW order. Away from half filling nesting no longer 
exists, however we still expect strong (tt, tt) antiferro- 
magnetic fluctuations. Using the arguments presented 
in Sec. II and Sec. Ill we expect SDW driven supercon- 
ducting pairing satisfying AkAk-i-q < 0. Both extended 
s-wave coskx + cos ky and d^^-y^-wave coskx — cos ky 
form factors satisfy this requirement. In the limit of 
small tz the nodes of cos k^ + cos ky lies close to both FS 
(Fig. (16)(a)). As the result dj.2^y2 will be favored. As 
tz increases, both the electron and hole pockets shrink. 
As they moves away from the node line of cos k^ + cos ky 
pairing (Fig. (18)(a)) the extended s-wave pairing will 
become favored over the d^^-y^ pairing because the lat- 
ter still produce nodes on the FS'*^. With the extended 
s-wave pairing the gap function on the two FS have oppo- 
site sign which is reminiscent of the gaps proposed in iron 
pnictides. We will follow the term used in the iron pnic- 
tide and call it s± SC. The above argument suggests that 
as a function of increasing tz there is a phase transition 
from c?-wave SC pairing to the s± pairing. A schematic 
phase diagram is shown in Fig. 15(b). 

In Fig. 16,17,18 we present the FRG result for this 
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model for tz = 0.5, 1.0 and 1.5. For tz = 0.5, as shown 
in Fig. 16(a), the FS intersects with the nodal line of 
cos kx + cos ky many time. In this shown in Fig. 

16(b), the Wi associated with d-wave pairing grows much 
faster than that of s-t-wave. The FRG form factor asso- 
ciated with these two types of pairing are shown in Fig. 
16(c) and (d), respectively. For tz — 1.0, the FS is shown 
in Fig. 17(a). Although the nodal line of cosfca, + cosfcj, 
no longer intersects the FS, a large portion of the hole 
FS lies close to the nodal line, hence will have a small 
gap. On the other hand, the nodal line of cos k^ — cos ky 
intersect with both two FS four times, giving rise to four 
nodes on each. It will be hard to tell which one will be 
favored by the kinetic energy. In fact, as shown in Fig. 
17(b), the Wi of both types of pairing are nearly degen- 
erate. As tz further increases, e.g., tz = 1.5 as shown in 
Fig. 18, the Wi associated with s± pairing dominates over 
all other pairing form factors. In short, our FRG calcula- 
tion confirms the earlier heuristic argument suggesting a 
transition between the dx2-y2 pairing and s± pairing. It 
would be very interesting to fabricate a double material 
with tunable interlayer hopping strength. For example 
it might be possible to realize this type of model using 
fermion optical lattices when the technical problem of 
cooling is overcome. 



V. THE RENORMALIZATION GROUP STUDY 
OF A LADDER VERSION OF THE IRON 
PNICTIDES 

In the cuprate physics, the study of the ladder ver- 
sion of the Hubbard model has contributed significantly 
to our understanding. Due to the quasi-one-dimensional 
nature of the system the Hubbard ladder is amenable to 
more exact treatment. Examples include density matrix 
renormalization group (DMRG), perturbative RG plus 
bosonization etc. Motivated by the above we perform 
the perturbative RG study for the two-leg ladder version 
of Eq. 5. The formulation, terminology and notations 
used in this section closely follows those of Ref. p^]. The 
readers can consult it for more details. 

The greatest simplification of the ladder RG, as com- 
pared to the FRG, is the fact that FS become Fermi 
points. As a result the functional RG equation become 
a group of coupling constant recursion relations. These 
equations can be integrated numerically to yield the cou- 
pling constants as a function of the energy cutoff A. As 
we shall show in the following, the qualitative nature of 
the ladder RG results agree with the two-dimensional 
FRG ones. In particular we also get the spin singlet s± 
pairing. In addition, the same pairing mechanism ap- 
plies. 

In the following we consider a 2-leg ladder running in 
the X direction. The y direction is consisted of two chains, 
and periodic boundary condition is applied in both direc- 
tions. The Hamiltonian we use is the same as Eq. 5. The 
only difference is the first Brillouin zone is changed from 
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FIG. 19: (color on-line) (a) The FS of iron pnictide and the 
five pairs of Fermi points for the iron pnictide ladder, (b) The 
Wi{A) associated with all singlet SC pairing. The red line is 
the Wi associated with the leading SC instability. 



-TT < kx.y < n to 



'TT < kx < TT, ky = 0, TT. Due to the 



discrete value of ky the five FS becomes five pair of Fermi 
points (Fig. (19)). Fortunately most of the high symme- 
try points on the 2D FS are captured. Following Ref.p^] 
we label these Fermi points as ni?, nL where n = 1, . . . , 5 
and R, L denote right and left movers respectively. It is 
worth to note that for the electron FS R, L correspond 
to the right and left members of the Fermi pair, and 
for the hole FS, R,L denotes the left and right mem- 
bers of the Fermi pair. The reversal is due to the fact 
that the dispersion around the hole FS is the opposite 
of that around the electron ones. More specifically we 
have li?,lL = (ipfc^i,0), 2i?, 2L = (ipfcF2,0), 3i?, 3L = 
(TfcF3,-7r), 4i?,4L = (±fc_F4,-7r), 5R,5L = (7r±fcF5,0), 
where 1,2,3 label the hole FS while 4,5 label the elec- 
tron FS. Associated with each Fermi pair there are the 
"slowly varying" chiral fields ijjRna and ipLna where a 
is spin index. The kinetic part of the Hamiltonian, lin- 
earized around the Fermi points, is given by 



5 



where vpn is the Fermi velocity associated with the n-th 
Fermi pair. 

In the following we assume there is no nesting (hence 
no Umklapp process exists). Under this condition, the 
set of momentum conserving four-point vertices have the 
same form as those in a 5-leg one band Hubbard ladder 
at generic dopings"^^. These interactions are given by 

TT. _ i~p jR jL _ ~a jR , jL \ 

where 



a a.l3 



Rmfi 



L 711/3 



Q,/3 
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with a being the three Pauh matrices. 

Since symmetry requires S^'^ = c^^ and j^^ = /^'^, 
we have a total of 2 x (15 + 10) = 50 coupling constants. 
In addition the ladder RG starts from a very small cutoff, 
so there is no "projection error" as in the FRG^®'^^. 

The RG equations can be easily obtained from Eq. 
(3.7)-(3.10) of RefP] 



In the above 



nra ^nra 



T^ivpn + VFm) ' 



(40) 



tp 

J mi 



J nra \JnmJ ' ^nm^nm nmJ 

f.P — o(rP fP A- —r'^ ~1 



k 

\ ^nmJnm ' ^nmJnm nrnJnmJ 

~ (^nm,k {c^nk^km + '^nk'^km + l^'^nk'^km) 



and 



^nm,k 



{Vpn + VFk){vFk + VPrn) 
2vFk{vFn + VFm) 



(41) 



where m n, and 

cP = (c" )^ + —(d^ )'^) 



k 

— 'y^,'^nn,k{c^k'^kn + ^nk^kn + l^^nk^kn)- 



Moreover the dot in e.g., denotes the derivative with 
respect to the logarithm of the energy cutoff. 

The differences between the 2-leg pnictidc ladder the 
the 5-leg one band Hubbard ladder lie in the values of the 
Fermi velocities vpn and the initial values of the coupling 
constants. Specifically the initial coupling constants are 
determined by the parameters in Eq. 5 and the band 
wavefunctions as follows 
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In the above equation a,b = 1, . . . , 5 label the five d- orbitals of Fe, UR^^naj a= 1, ...5 are the five orbital com- 
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FIG. 20: (color on-line) Left panels: the schematic represen- 
tation of the strongest scattering processes (their scattering 
amplitude is at least two orders of magnitudes larger than all 
other scattering amplitudes) for the 2-leg pnictide ladder at 
low energies: (a) C24 (solid line) and C22 (dashed line) (b) C24 
and C22 (dashed line) (c) /24 (solid line). Right panels: the 
RG flow associated with the processes associated with the 
right panels. 



ponents of the Bloch wave functions associated with the 
R/L members of the nth Fermi pair. Finally (a <-> b) 
denotes the exchange the subscript a and b in the term 
right before it. 

We stress here that most of the above coupling 
constants have the dual character described in ear- 
lier sections. For example, the term £24 J^Jik 



T,a.,t3^2i^R2o.'^[2f}'^ Lipoma. Can be viewed as both 
Cooper scattering and SDW scattering between hole 
pocket 2 and electron pocket 4. In this sense, the lad- 
der RG overemphasis these special channels. As a result 
the ladder RG retains the key ingredient which allow the 
SDW scattering to drive other ordering tendencies. 

Numerically we found that all the coupling constants 
diverge at a finite cutoff. Moreover the coupling con- 
stant associated with two Fermi pairs, one from the elec- 
tron FS and the other from the hole FS, always domi- 
nate over others. For example, with the undoped band 
structure and initial couplings determined by the same 



set of parameters, Ui = 4, C/2 = 2, Jh = 0.7, used in 
Sec. II we found that the coupling constants associated 
with Fermi pair 2 (from the large hole pocket around 
(0,0)) and 4 (from the electron pocket around (0, — tt)) 
are larger than others by more than two orders of magni- 
tude. More specifically the coupling constants with the 
largest magnitudes are 



-24' 



-22' 



^44' 



^44 ' J 2' 



24- 



(42) 



The schematic representation of these scattering pro- 
cesses and their RG flow are plotted as Fig. 20. 

As mentioned before C24 has the dual character of SDW 
and Cooper scattering (solid line of Fig. 20(a)), it grows 
to large positive value hence drives SC pairing with re- 
versed sign between electron and hole pocket. C24 has the 
dual character of Cooper scattering and CDW (solid line 
of Fig. 20(b)), it also flows to large positive value, and 
contributes to both s± SC and imaginary CDW order. 
Here we note a little departure from the 2D physics in the 
sense that in 2D it is the dual SDW and CDW scattering 
that drives imaginary CDW. Here due to the omission of 
the Umklapp process LL — > RR between n = 2 and n = 4 
the process in Fig. 7(a) is disallowed. C22 and £44 are intra- 
pocket Cooper scattering (dashed line in Fig. 20(a,b)), 
they are driven by the inter-pocket Cooper scattering 
(which favors the s± pairing) and become large negative 
at the lowest energies. /24 is the shared channel between 
SDW and forward scattering (solid line in Fig. 20(c)), it 
flows to large positive and causes PI as discussed in Sec 
II. It is interesting to observe that all important compet- 
ing orders found in 2D FRG are reproduced in this ID 
ladder RG. 

Pairing form factors are studied in the same way as in 
Sec. II. In this case Eq. (15) are 5x5 matrices (n,™ = 
1,...,5) 



Vsc,t(?^,"^, A) 



4 ™ 



and negative eigenvalues Wi indicate pairing instability. 
We found that the triplet channels are much weaker than 
the singlet ones, and among the five singlet channels only 
one of them become negative upon RG flow Fig. 19(b) 
Eigenvector(form factor) corresponding to this channel 
has opposite signs between the hole pocket 2 and elec- 
tron pocket 4, consistent with the s± pairing found in 
2D FRG, but has nearly vanishing weights on the other 
Fermi points. 



VI. CONCLUSION 

Motivated by the high Tc superconductivity of iron 
pnictide, we study the pure electronic pairing mecha- 
nism in both the cuprates and the pnictides from an 
itinerant point of view, where the FRG method becomes 
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very useful. The main message of this paper is a view 
point we obtained from the FRG studies, namely, in a 
system with strong antiferromagnctic correlation such as 
the cuprate and the iron pnictide, there are generically 
SC, PI and orbital-current (CDW) ordering tendencies. 
These ordering tendencies are triggered by a set of dual- 
character scattering processes associated with SDW-SC, 
SDW-PI and SDW-CDW, From these dual processes, 
other generic SC, PI and CDW scattering grow. Thus 
it is through these dual processes that SDW fluctuation 
drives all the other instabilities. Since these SDW related 
dual processes all have positive amplitudes, it imposes 
strong constraint on the symmetry of the form factor as- 
sociated with the induced order parameters. The details 
as summarized in the table below. 

Other significant results presented in this paper are 
also summarized as follows, (i) The origin of the SC 
gap anisotropy in iron pnictides. We argue the SC gap 
anisotropy in the iron pnictide is due to the variation of 
the orbital content of the Bloch wave function along the 
FS. (ii) We show that the Pomerachuk instability in the 
iron pnictide tends to shrink both the electron and hole 
pockets. It can give rise to a change of the FS topol- 
ogy of the electron pockets, (iii) We demonstrate the 
validity of the Ji — J2 model in describing the lowest 



energy effective Hamiltonian for the iron pnictides, and 
determine the range of Ji — J2 for which this model is 
valid. Wc commented on the inadequacy of the J\ — Ji 
model in the intermediate energy regime where compet- 
ing order is manifested, (iv) We discuss the ordering 
tendencies of the cupratcs including charge stripes under 
the same framework, (v) Wc show that the double-layer 
Hubbard model nicely interpolates between the physics 
of the cupratcs and the physics of the pnictides as a func- 
tion of the inter layer hopping parameter. We predict the 
existence of a quantum phase transition where the super- 
conducting pairing changes from (ij.2_,y2-wavc to the s± 
symmetry, (vi) We present the rcnormalization group re- 
sults for a 2-leg ladder version of the iron pnictides Hamil- 
tonian. We show that the results are qualitatively consis- 
tent with those obtained from two-dimensional functional 
renormalization group calculations. 
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TABLE I: Summary of spin density wave fluctuation driven superconducting pairing, Pomeranchuk instability and charge 
current density wave 
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